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ABSTRACT 

Two methods are used to observationally determine the solar radius: One is the observation 
of the intensity profile at the limb, the other one uses f-mode frequencies to derive a 'seismic' 
solar radius which is then corrected to optical depth unity. The two methods are inconsistent 
and lead to a difference in the solar radius of ~0.3Mm. Because of the geometrical extention 
of the solar photosphere and the increased path lengths of tangential rays the Sun appears to 
be larger to an observer who measures the extent of the solar disk. Based on radiative transfer 
calculations we show that this discrepancy can be explained by the difference between the height 
at disk center where T5000 = 1 (tross — 2/3) and the inflection point of the intensity profile on the 
limb. We calculate the intensity profile of the limb for the MDI continuum and the continuum at 
5000 A for two atmosphere structures and compare the position of the inflection points with the 
radius at T5000 = 1 (tross = 2/3). The calculated difference between the 'seismic' radius and the 
inflection point is 0.347 ± 0.06 Mm with respect to T5000 = 1 and 0.333 ± 0.08 Mm with respect 
to tross = 2/3. We conclude that the standard solar radius in evolutionary models has to be 
lowered by 0.333 ± 0.08 Mm and is 695.66 Mm. Furthermore, this correction reconciles inflection 
point measurements and the 'seismic' radii within the uncertainty. 

Subject headings: astrometry — Sun: fundamental parameters — Sun: photosphere — radiative transfer 



1. Introduction 

The solar radius is a key parameter for under- 
standing the Sun's evolution and many astrophys- 
ical applications. Being our closest star, the Sun 
can be considered as our laboratory, allowing us 
to determine its radius with high precision. Over 
the past 25 years various groups were dedicated to 
determine the solar radius by applying basically 
two approaches. First, there is the direct measure- 
ment of the photospheric radius. In earlier times 
its visual determination was carried out by an ob- 
server. Later, employing CCD-imaging devices, 



the visual radius determination was replaced by 
deriving the inflection point to the intensity profile 
of the limb. The error of both methods has been 
investigated in detail bv lLaclare et al. I (|l996h who 
derive a solar radius of 959" .60, which agre es well 
with 959" .63 ± 0".10 or 695.99 ±0.07 Mm (|Allenl 



19761 ). commonl y used for calibrating solar evolu 



tiona r y models dBrown &: Christensen-Dalsgaar' 
19981: iTurck-Chieze fc Murdinl l2000h . Second, 



Schou et all (|l997t k lAntial (|l998h determine the 



helioseismic radius by means of f-mode frequen- 
cies o btained from the M ichelson Doppler Im- 
ager (|Scherrer et al.l Il995l MDI) on board the 
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ES A/NASA spacecraft SOHO. Due to the high 
quality of the MDI observations the helioseismic 
determination is very precise and leads to a ra- 
dius that is approxim ately . 3 Mm smaller than 
the st andard value by All enl (|l976l ). ISchou et al.l 
(|l997l ) conclude that the standard solar radius 
has to be decreased by this va lue. Finally, 



Brown fc Christensen-Dalsgaardl ( 1998h deter- 
mine the solar radius from observations obtained 
with the High Altitude Observato ry's Solar Diam- 
eter Monitor (jBrown et al.lll982l HAO SDM) by 
applying a model of the data reduction procedure. 
The authors point out that there are significant 
uncertainties associated with the currently used 
radius value which are related to the definition of 
the solar limb adopted in the radius determina- 
tions and the reduction of the measured value to 
the photosphere. The authors indicate a height of 
about 0.3 Mm which 'perhaps' explains the radius 
correction inferred from the f-mode frequencies. 
In the present paper we determine the difference 
between the height at r 50 oo = 1 at disk center 
and the distance of the inflection point from ra- 
diative transfer calculations for two wavelengths 
and atmosphere structures. 

2. Calculations 

With the radiative transfer code COSI we cal- 
culate the solar limb function for the MDI contin- 
uum intensity and the continuum at 5000 A. COSI 
is a combination of a model atmo sphere code in 



spheri cal symmetry, developed by iSchmutz et al 



(1989) and the spe ctrum synthesi s program SYN- 
SPEC, goin g back to Hubenv (1981) and further de- 
veloped by lHubenv & Land (|l992al |bl). The adap- 
ti on to the solar atmosph e re ha s been published 
Haberreiter fc Schmutd ( 20031 ). The model at- 



mosphere code calculates the NLTE population 
numbers for a set of specified atomic levels by 
solving the radiative transport equations simulta- 
neously with the equations for statistical equilib- 
rium. The radiative transfer is solved along rays 
parallel to the central ray incident on a spherical 
distribution of the p hysical atmosphere structure 
( Mihalas et al.fl975h . A spherical geometry allows 
to calculate geometrically correct the emerging in- 
tensity at the limb and line of sights beyond the 
solar limb. 

We base the model calculations on two model 



atmosphere structu res for the q uiet Sun, the 
LTE structure by ( Kurucd 1991 . hereafter Ku- 
rucz91) and the NL TE structure Model C by 
Fontenla etaH (|2007l hereafter Fontenla07). The 
radius at disk center where T5000 = 1, is set to the 
solar radius derived from s eismic obser v ations , 
i? =695.68 Mm as given bv ISchou et~all (|l997h . 
The high resolution Ni line profile has been con- 
volved with the five MDI filter functions F0 to F4 
used to determine the actual continuum intensity 
observed by MDI. 

3. Results 

Fig. Q] panel (a) shows the mean intensity varia- 
tion for the MDI continuum at the solar limb ap- 
plying all five MDI filter functions and the contin- 
uum at 5000 A calculated with the LTE and NLTE 
model atmosphere structures. Distance zero refers 
to the 'seismic' solar radius i?0=695.68 Mm, i.e. 
T 5000 = 1 at disk center. For the exact determina- 
tion of the inflection point we calculate the deriva- 
tive dl/dx (x being the outward coordinate) of all 
four limb profiles, shown in Fig. [T] panel (b). The 
two atmosphere structures agree well in the pho- 
tosphere and lead to very similar intensity profiles. 
The derivative of the intensity is always negative 
and shows a local minimum at the inflection point. 
The distances of the local minima for all four limb 
profiles are given in Table[TJ The mean value of 
the distance for the two atmosphere structures for 
A = 5000 A is 0.347±0.06 Mm, where we used the 
difference between the two atmosphere structures 
as criterion to derive the error introduced by the 
atmosphere structures. This re sult is close to the 
value of 0.370 Mm estimated bv lWittmaiml ll 19741) 



based on the atmosphere structure bv lllolwege] 
(1 1967ft 



In evolutionary models the radius is typically 
defined as the point where the temperature equals 
the effective temperature, i.e. where tr oss = 2/3. 
From our calculations we find that the height 
where tr oss = 2/3 is in the mean 0.014 Mm fur- 
ther out than where T5000 = 1- Thus the cor- 
rection of an inflection point measurement to the 
'seismic' radius is 0.333 Mm. This correction ex- 
plains the discrepancy between the modeled and 
observed f-mode frequencies. We conclude that 
the standard radius used in the evolutionary mod- 
els has to be corrected by value this value, leading 
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Fig. 1. — Panel (a) gives the intensity profile calculated with the Fontenla07 and Kurucz91 solar atmosphere 
structures for the 5000 A continuum (dotted and dashed-dotted) and the MDI continuum (solid and dashed) . 
Panel (b) shows the derivative of the four intensity profiles. The local minimum is the position of the inflection 
point, which is 365 km for MDI and 347 km for 5000 A (see TableQ]). Note that the distance to tr, oss = 2/3 
is 14 km less. 
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to 695.66 Mm. 

Table 1: Distance z (Mm) between the inflection 
point and the r 50 oo = 1 radius (i? Q =695.68 Mm) 
for two different atmosphere structures and wave- 
lengths. The distance to tr oss = 2/3 is in the 
mean 0.014 Mm less. 





T5000 — 1 


Atmosph. struct. 


MDI cont. 


5000 A cont. 


Kurucz91 


0.367 


0.344 


Fontenla07 


0.370 


0.350 


Az 


0.003 


0.006 


z 


0.368 


0.347 



Table 3: Correction of 0.333 Mm applied to the 
mean of the inflection point measurements. The 
corrected radius value is compared with the mean 
of the 'seismic' radii given in Table [2) referring to 
tross = 2/3. The result is consistent within the 
uncertainty. 

Radius (Mm) 

Rip 695.901±0.098 

Correction _ 0.333±0.008 

Corrected R Q 695.568±0.098 

R Q 695.658±0.140 

AR 0.090±0.171 



Discussions 



For consistency we investigate whether recent 
inflection point measurements are compatible with 
the mean of 'seismic' radii relating to 7r oss = 2/3, 
i.e. where the temperature equals the effective 
temperature, by applying this correction term. 
Table [5] gives the solar radii based on ii) IP mea- 
surements, ii) 'seismic' radii. The mean values 
deviate by ~0.24 Mm. In Table[3]we correct the 
mean radius values of the inflection point obser- 
vation by 0.333 Mm and compare the result with 
the mean 'seismic' radius, leading to a difference 
of 0.09±0.17 Mm. We conclude that the inflec- 
tion point measurements and the 'seismic' radius 
can be reconciled by applying the correction of 
0.333 Mm. 

Our calculations indicate that the position of 
MDI continuum inflection point is by 0.02 Mm 
further out than for 5000 A. T he radius value de- 
termined bv lKuhn et al. I (|2004l ) (see Table© refers 
to the inflection point of the intensity profile ob- 
served by MDI. The difference of 0.02 Mm can- 
not explain the so mewhat lower radius value by 
Kuhn et all lj2004h . 

We emphasis that we do not consi der any varia- 
tion of the solar diame t er wit h time dEmilio et all 
20071: iThtiillier et all [200l iKuhn et all |2004 
Toulmondd Il997t IParkinson et all Il980h . which 
however could explain some of the deviation of 
the radius values. Furthermore, our considera- 
tions are with respect to the quiet Sun. 



Conclusions 



From radiative transfer calculations with the 
spherical code COSI we determine a distance of 
0.333±0.008 Mm (0.347 ± 0.06 Mm) between the 
height where tr ss = 2/3 (T5000 = 1) at disk cen- 
ter and the position of the inflection point of the 
5000 A intensity profile at the limb. This correc- 
tion explains the differences between the f-mode 
frequencies derived from models calculations and 
observations. Therefore, we conclude that stan- 
dard solar radius currently used in evolutionary 
models has to be corrected by 0.333 Mm and is 
695.66 Mm. 

Furthermore, the correction of 0.333 Mm rec- 
onciles the radius values derived from inflection 
point measurements and the 'seismic' radius by 
0.09±0.17 Mm. Due to the increasing precision 
of the instruments it becomes necessary that the 
community agrees on a generally binding defini- 
tion of the solar radius and that values derived 
via different techniques are to be corrected ac- 
cordingly. We propose that the solar radius refers 
to i?© = i?(TR,oss = 2/3). Thus, we suggest 
that the value of the solar radius determined from 
inflexion point measurements has to be lowered 
by 0.333 Mm to be comparable with 'seismic' 
radii. A detailed study of the position of the 
wavelength and solar cycle dependent inflection 
point will be crucial for upcoming observations 
with the SODISM instrument on board PICARD 



(|Thuilher et alj 120061 ). 
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Table 2: Solar radii and their mean values Rq derived by means of i) the standard solar radius and recent 
i nflect ion point measurements ii) values referring to the 'seismic' radius. As the result given by IWittmann 
(Il997h more than 8a different from the mean we excluded this value from our considerations. 



Authors 


R (Mm) 


R ( 


arcsec) 




Description 


R & (Mm) 


i) IP measurements 














Allen (19761 


695.990±0.070 


959" 


63±0" 


02 


standard value 




Sofia et al. (19941 


695.917±0.043 


959" 


53±0" 


06 


SDS 




Neckel 19951 


695.982±0.022 


959" 


62±0" 


03 


McMath ST 




Laclare et al. (19961 


695.830±0.007 


959" 


41±0" 


01 


OCA astrolab 




Kuhn et al. (20041 


695.740±0.110 


959" 


29±0" 


15 


MDI 




Emilio & Leister (20051 


695.946±0.029 


959" 


57±0" 


04 


SP astrolab 


695.901±0.098 


Wittmann ( 19971 


696.715±0.029 


960" 


63±0" 


04 


Obs. del Teide 




ii) 'seismic ' radii 














Schou et al. (19971 


695.680±0.030 


959" 


20±0" 


01 


MDI f-modes 




Antia (19981 


695.787± 


959" 


35± - 




GONG f-modes 




Brown & CD (1998) 


695.508±0.026 


958" 


97±0" 


04 


HAO SDM 


695.658±0.140 



the Project No. 200020-109420. 
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